o 

(N 

< 



First-principles study of heat transport properties of 

graphene nanoribbons 

Zhen Wah Tan,"'''* Jian-Sheng Wang,*)^ and Chee Kwan Gan*'* 

Department of Physics, Cornell University, 109 Clark Hall, Ithaca, New York 14853-2501, 
U.S.A., Institute of High Performance Computing, 1 Fusionopolis Way, #16-16 Connexis, 



E-mail: 'ganck@ihpc.a-star.edu.sg 



Abstract 



a Singapore 138632, Singapore, and Center for Computational Science and Engineering, and 

czi Department of Physics, National University of Singapore, Singapore 117542 

e 

o 

We use density-functional theory and the nonequiUbrium Green's function method, as well 

(N 

K^ as phonon dispersion calculations to study the thermal conductance of graphene nanoribbons 

^ with armchair and zigzag edges, with and without hydrogen passivation. We find that low- 

frequency phonon bands of the zigzag ribbons are more dispersive than those of the armchair 

^ ribbons, and that this difference accounts for the anisotropy in the thermal conductance of 

. . graphene nanoribbons. Comparing our results with data on large-area graphene, edge effects 

'k> are shown to contribute to thermal conductance, enhance the anisotropy in thermal conductance 

c3 of graphene nanoribbons, and increase thermal conductance per unit width. The edges with 

and without hydrogen-passivation modify the atomic structure and ultimately influence the 

phonon thermal transport differently for the two ribbon types. 
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Introduction 

Graphene nanoribbons (GNR), or strips of planar graphenejl^ exhibit several novel properties such 
as large magnetoresistance'^ and high electron mobility's which prove to be useful in nanoscale 
electronics. Graphene is known to have extremely high thermal conductivityl^ that depends on the 

flake size's] and number of atomic planes.^ GNR are expected to retain good thermal properties. 

pn 
This, together with the availability of simpler fabrication techniques than carbon nanotubes, has 

caused GNR to gain much attention in recent years due to their potential applications in modern 
devices. Experimental and theoretical work have concluded that GNR can be metallic or semi- 
conducting,'^]^ and the energy gap varies with ribbon width'^mnH ^^^ edge orientation.-^ These 
properties make semiconducting GNR an attractive alternative channel material'^ capable of pro- 
ducing smaller devices than those achievable with silicon. Metallic GNR can also be used as 
interconnects, ^^ and having circuits made entirely of GNR holds the possibility of reducing or 
eliminating contact resistance. ^^ For these reasons, various methods for fabricating GNR have 
been developed,' ^*^ ' ^^ * ^^ ' including scanning tunneling microscope lithography that offers precise 
control of the structure of ribbons. Recent developments have also enabled the production of 
sub-10-nanometer GNR,'^and greater control over edge geometries.!^ Prototype GNR field-effect 
transistors have also been developed and characterized.'^^'^Il 

While the electronic properties of GNR have been explored in depth, thermal transport in 
GNR — an important consideration for thermoelectric performance and thermal management — has 
also gained more attention recently .'22n23 Several studies have been made on the phonon disper- 
sions'22n20| Qf gnr^ while more recent articles have focused on phonon-mediated thermal trans- 
portpH of GNR. For example, Sevin9li and Cuniberti have investigated the thermoelectric perfor- 
mance of edge-disordered GNR.'^^ Hu et al.^ observed that ZGNR have higher thermal conduc- 



tance than AGNR of comparable widths, and attributed this to different phonon scattering rates 
intrinsic to the ribbon geometry. Several other studies have also confirmed the anisotropy in ther- 
mal conductance of GNR,'^^'^! ^^^ ^ similar anisotropy has been predicted by Jiang et al.^ in 
graphene. In this Letter, we employ density-functional theory (DFT) that accurately deals with the 
atomic and electronic structures, the nonequilibrium Green function method that has been shown 
to be useful in studying the electron ^^ and heat transport through nanostructures,'23 and also the 
standard theory of lattice dynamics to account for phonon transport through GNR. We investigate 
how the atomic structures of GNR affect the phonon dispersions and ultimately thermal transport, 
and explain the anisotropy in thermal conductance of GNR. 

Methodology 
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Figure 1: Hydrogen-passivated graphene nanoribbons with (a) armchair edges (AGNR-A^), and (b) 
zigzag edges (ZGNR-A^). A primitive unit cell is marked out in each case, and the ribbons are 
periodic in the a-direction. The number A^ denotes the number of C atoms in the Z^-direction, and 
the ribbon width W is given by the maximal distance between C atoms in the Z?-direction. 

Because systems of sub-nanometer widths will become important with the miniaturization of 
GNR-based devices, we study the thermal conductance in narrow, pristine GNR by combining 
first-principles density-functional calculations and the ballistic nonequilibrium Green's function 
(NEGF) method.'22Ell j^ this work, we investigate the thermal conductance of armchair GNR 



(AGNR-//) of widths // = 3,4,5,6 (refer to [figure][l][]l i for the meaning of A^), as well as 



zigzag GNR (ZGNR-A^) of widths A^ = 2,3,4,5 ( |[figure][l][]"T» . We obtain the force-constant 



matrix via DFT calculations implemented in SIESTA.'^ We use periodic boundary conditions on 
the orthorhombic supercell consisting of nine primitive unit cells: The ribbons are periodic in the 



a-direction (see Figure [figure][l][]l ), and each ribbon is separated from its nearest neighbors by 

o 

16 A of vacuum in the b- and c-directions. To facilitate NEGF calculations, we partition the su- 
percell into three blocks with three primitive unit cells per block. These three blocks are for the 
central junction, the left and right lead regions as described in Ref.'^Sl We note that the same force 
constants required for NEGF calculations can be readily used to calculate the phonon dispersion 
curves. 

For DFT calculations, the local-density approximation with the exchange-correlation func- 
tional due to Perdew and Zunger is used. We use a rather fine mesh cutoff of 400 Ry. As 
demonstrated by Son et al.'^ and Gan et a/. ,33 spin-polarization effects are particularly important 
in ZGNR, therefore we use spin-polarized calculations for ZGNR in this work. Nonspin-polarized 
calculations are used for AGNR. 

We first perform ionic relaxation of GNR using a conjugate-gradient method. To obtain the best 
possible relaxed structure while ensuring that the conjugate-gradient minimization converges, we 
use a force tolerance of 0.001 eV/A. To construct a force-constant matrix, we sequentially displace 
each atom from its equilibrium position in the a-, b- and c-directions by a distance of 0.015 A, 
and evaluate the forces acting on all atoms as a result of each displacement.'^ A central finite 
difference scheme is used to evaluate the force-constant matrix, which is required by the NEGF 
and the phonon calculations. 

Using the ballistic NEGF method,'22ESl ^ve then calculate the phonon transmission coefficient 
f for each GNR system. Thermal conductance is evaluated using the Landauer formula: 

<T(r) = /^°°gf.cof(co)|^, (1) 

where the occupation distribution function f{co,T) = l/(e^'^/^B^ — 1). It is important to note 
that ?? implies that at low temperatures, low-frequency modes are the dominant factor in thermal 



transport since the derivative df/dT diminishes rapidly with increasing (O. 



We find the thermal conductance of GNR to be independent of ribbon length in the ballistic 
regime, as noted previously by Lan et alF^ This is consistent with the theory of ballistic transport 
based on the Landauer formula, where diffusive behaviors cannot be addressed. The ballistic as- 
sumption is valid for small systems, such as the ones studied in this work, where the system size is 
much smaller than the graphene phonon mean free path of ~ 775 nm at room temperature's How- 
ever, in the case of large-area graphene, diffusive (Umklapp-limited) scattering plays a significant 
role in reducing thermal conductivity, as determined theoretically^^ and experimentally.S^IBl 

Results 
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Figure 2: (color online) Thermal conductance o (at 300 K) of GNR as a function of ribbon width 
W, for ZGNR-A^ with A^ = 2, 3, 4, 5 and AGNR-A^ with A^ = 3, 4, 5, 6. Data shown are for ZGNR 
without H-passivation, ZGNR with H-passivation; AGNR without H-passivation, and AGNR with 
H-passivation — linear fits have correlation coefficients 0.998, 1.000, 0.882 and 0.978, respectively. 
Gradients of the linear fits are 1.276, 1.203, 1.119 and 0.771 W/(m K), while vertical intercepts are 
0.4523, 0.4906, 0.178 and 0.3967 nW/K, respectively. 



[figure] [2] []2 shows the thermal conductance (at 300 K) of various GNR as a function of rib- 



bon width W . In the limit of large width, the conductance of a ribbon is proportional to its width, 
and our data exhibit a similar trend: For ZGNR with and without H-passivation, and for AGNR 
with H-passivation, the thermal conductance increases linearly with width, with high correlation 
coefficients of 0.978 and above. While thermal conductance of AGNR without H-passivation gen- 
erally increases with temperature, it does not follow a strictly linear trend. On the other hand, on 



extrapolating the linear trends toW = A, we see positive contributions to the thermal conduc- 
tance (~0.5 nW/K) that stem from edge effects. This is consistent with the findings by Xu et al.p^ 
which show a sharp increase in thermal conductance per unit width as W decreases from 1 nm to 



0.5 nm, where edge effects becomes more evident, [figure] [2] [] 2 also shows that ZGNR generally 
have higher conductance than AGNR of comparable widths. The main reason for this is due to the 
larger transmission coefficients f{v) for the ZGNR compared to AGNR of comparable widths, as 



shown in [figure][3][]3 This observation is explained in the following paragraphs. 
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Figure 3: (color online) Comparison between phonon dispersions (left panel) and transmission 
coefficients (right panel) for AGNR-5 and ZGNR-3, (a) without and (b) with H-passivation. 
Imaginary-frequency modes are indicated by v < cm^^. The two ribbons have comparable 
widths: W^agnr-s = 4.942 A, Wzgnr-3 = 4.993 A. 

According to the NEGF formalism for ID systems,S2MII ^^g transmission coefficient T{v) 



is the number of phonon modes at any frequency v = (o/2n. This can be readily seen in [fig 



ure] [3] []3 where we compare the phonon dispersions with the transmission coefficient for AGNR- 
5 and ZGNR-3, with and without H-passivation. We see that T{v) essentially equals to the number 
of bands at frequency v, except at very low frequencies. We expect that the discrepancy at low 
frequencies between f (v) from NEGF calculations and that deduced from phonon bands is caused 



by intricate differences in the numerical treatments for these two methods. For example, in our 
NEGF calculation, it is assumed that two atoms cease to interact with one another (i.e., the force 
constant between them is set to zero) once they are are separated by more than three primitive 
cells. However, we do not truncate force-constant matrix elements in the phonon calculations. The 
slightly imaginary modes near F from the phonon calculations suggest that the GNR is unstable 
with respect to long wavelength periodic distortion, a fact that is consistent with the compres- 
sive edge stresses that are inherent in the GNR edges. ^'^ Despite small numerical uncertainties for 
extremely low frequency modes (which may be eliminated by adopting a larger supercell at the 
expense of incurring larger computational costs), the data obtained for f{v) is sufficiently accu- 
rate for conductance calculations. Since f{v) is equal to the number of phonon modes present, 
dispersive phonon bands would generally lead to larger values of T{v) than less dispersive bands. 



The dispersion relations in [figure] [3] [] 3 indicate that low-frequency phonon bands in ZGNR- 
3 tend to be more dispersive than those in AGNR-5. The ribbons have comparable widths, with 
V7zGNR-3 = 4.993 A, and Wagnr-5 = 4.942 A. The dispersive bands in ZGNR-3 would give rise to 
higher f{v) in the low-frequency regime: low-frequency modes play a dominant role in thermal 
conductance at room temperature. Thus, the fact that low-frequency bands in ZGNR-3 are more 
dispersive than those in AGNR-5 explains the anisotropy in thermal conductance, even in the 
ballistic limit. 

Since the transmission coefficients are intimately related to the phonon dispersions, we have 
systematically studied the phonon dispersion of the ribbons of different width and orientations 
(i.e., AGNR-A'^ and ZGNR-A^), with and without hydrogen passivation, where the results are sum- 



marized in [figure] [4] []4 The dispersion trends, cut-off frequencies, and C-H stretching mode 
frequencies agree very well with previously published data calculated using REBOII and MO/8 
simulations.'^^'^oi Also, DFT calculations by Gillen et a/. I^SI yielded very similar phonon dispersion 
trends for ZGNR, with cut-off frequencies of ~1650 cm^^, and C-H stretching mode frequencies 
of ~ 3100 cm^^, while we obtain ~1700 cm^^ and ~3120 cm^^ respectively. 

It is interesting to note that in most of the test cases, the lowest phonon modes are found to 
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Figure 4: Phonon dispersion plots for AGNR-4, 5, 6 and ZGNR-4, 5 (a) without H-passivation and 
(b) with H-passivation. Imaginary frequencies are indicated by v < cm^^ 



soften and have imaginary frequencies — an indication of instability — near T. These imaginary- 



frequency modes (denoted by v < cm^) are shown in [figure] [4] []4 The lowest modes are edge 
waves with out-of-plane displacements, and the unstable edge rippling generally occurs over a 
range of about 10% the Brillouin zone size, corresponding to wavelengths larger than 2 nm (4 nm) 
for ZGNR (AGNR). This result is consistent with a study by Shenoy et al.j^ which indicates that 
GNR form edge ripples with wavelengths of about 8 nm due to edge stresses. 



The phonon plots in [figure] [4] []4 show that in general, ZGNR have low-frequency bands that 
are more dispersive than those in AGNR. The reasoning for ZGNR-3 having higher thermal con- 
ductance than AGNR-5 can thus be applied also to relate the anisotropic behavior of GNR: More 
dispersive low-frequency bands in ZGNR result in larger transmission coefficients f{v) as com- 
pared to AGNR in the low-frequency regime, thereby giving rise to higher thermal conductance in 
ZGNR. 

A similar anisotropy in thermal conductance in graphene sheets has been pointed out by Jiang 
et al. , -^^ using the valence force field model. The study found a thermal conductance per unit width, 
o /W = 0.348 and 0.352 W/(m K) in the armchair and zigzag directions, respectively. Similar find- 



ings have been obtained by Hu et alJ^ and Guo et alJ^ through molecular dynamics simulations 
for GNR with larger N. Hu et al. attributed the difference in thermal conductance between AGNR 
and ZGNR to the presence of different phonon scattering mechanisms. However, since we ob- 
serve that the anisotropy persists in the ballistic regime, we suggest that the anisotropy in thermal 
conductance could also be due to the anisotropy in phonon dispersion for GNR. 

On the other hand, the fitted gradients of our o versus W plots for H-passivated ribbons are 
0.771 and 1.203 W/(m K), respectively. This suggests that edge effects increase the o /W ratio, 
and magnify the anisotropy in thermal conductance for GNR. These trends indicate that arrays of 
narrow ZGNR may be most efficient in thermal management, whereas wide AGNR may optimize 
thermoelectric performance of graphene-based transistors, with its lower thermal conductance. 

Table 1 : Thermal conductance data of GNR. Thermal conductance at 300 K is a, in units of nW/K. 
The conductance shift due to H-passivation is A = (awithH — 0'w/oh)/<7w/oH- 



A^ 


AGNR 


ZGNR 




<^w/oH 


CwithH 


A 


<7w/oH 


<7withH 


A 


2 








0.797 


0.828 


3.9% 


3 


0.409 


0.570 


39.4% 


1.119 


1.096 


-2.1% 


4 


0.712 


0.719 


1% 


1.361 


1.358 


-0.2% 


5 


0.630 


0.757 


20.1% 


1.627 


1.599 


-1.7% 


6 


0.897 


0.875 


-2.5% 









In considering the effects of H-passivation on o, we define the thermal conductance shift by 



A = (awithH — 0'w/oh)/o'w/oH- Table [table][l][]l shows that while A stays below 4% for ZGNR, 



A can reach as high as 39% for AGNR. This disparity can be traced to the modification of phonon 
modes with H-passivation. 



[figure] [4] []4 shows that H-passivation induces little change to the dispersion relations for 



ZGNR, compared to AGNR. For ZGNR, despite the addition of six phonon bands (due to two 
extra H atoms), these bands have high frequencies (v > 700 cm^) and thus do not contribute sig- 
nificantly to thermal conductance at room temperature. On the other hand, for AGNR, we see that 
H-passivation causes more drastic changes to phonon dispersions. Besides having twelve more 
phonon bands due to the four H atoms, the phonon bands also generally shift towards lower fre- 
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quencies. To explain the latter, we note that H-passivation reduces the strong C-C triple bonds at 
the armrests of AGNRH2E0]jq ^ much weaker edge bond. However, in ZGNR, only one in two edge 
C atoms have dangling bonds, thus edge C-C bonds are weaker than double bonds, and relax only 
slightly upon H-passivation. Bond relaxation causes a reduction in force constants and phonon 
frequencies, and this change takes place more significantly in AGNR than in ZGNR. 

It should also be noted that for AGNR with and without H-passivation, the thermal conductance 
varies non-trivially with ribbon width, with larger variations observed in the cases without H- 
passivation. Previous studies on AGNR have shown that the variation of energy gapslSESEII and 
edge energies'^ with ribbon width can be grouped into three distinct families (with respect to A^ 
mod 3). However, we have been unable to verify such trends in thermal conductance for AGNR 
due to huge computational costs for systems of larger N. 

Conclusions 

We have studied phonon-mediated thermal conductance of graphene nanoribbons with armchair 
(AGNR) and zigzag (ZGNR) edges of sub-nanometer ribbon widths through density-functional 
calculations, a nonequilibrium Green's function method, and phonon calculations. ZGNR was 
found to have higher thermal conductance than AGNR of comparable widths, due to an anisotropy 
in the phonon dispersion for GNR: low-frequency bands in ZGNR are more dispersive than those 
in AGNR. Edge effects were found to produce a positive contribution (~0.5 nW/K) to thermal con- 
ductance, increase thermal conductance per unit width, and cause ZGNR to have a significantly 
higher thermal conductance per unit width (~1.2 W/(m K)) than AGNR (~0.77 W/(m K)). These 
facts suggest that narrow ZGNR can act as good thermal conductors for thermal management, 
while wide AGNR may be a better candidate for application in nanoscale devices. Thermal con- 
ductance of AGNR and ZGNR change differently with hydrogen-passivation at the edges: Edge 
C-C bonds in AGNR are relaxed by a greater extent than those in ZGNR, thus phonon frequencies 
in AGNR are more significantly reduced than those in ZGNR. It remains interesting to determine if 
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thermal conductance in AGNR follows the same trends as that for energy gap and edge energy, as 
this will yield further insight on thermal transport mechanisms in AGNR. This may shed some light 
in understanding the thermoelectric behavior of AGNR, and ultimately be useful for applications 
in nanoscale electronics. 

Finally, we note that the combination of density-functional theory, nonequilibrium Green's 
function method and phonon dispersion calculations reveal, for the first time, the peculiar heat 
transport properties that could be traced to the intricate interplay between the atomic structures, the 
transmission coefficients and the dispersiveness of the phonon bands. Our study provided a natural 
link between these quantities and explained the anisotropy observed in the thermal conductance of 
zigzag and armchair GNR. 
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